On the Reflection of Electrons by Metallic Crystals 
By L. A. Mac COLL 

This paper gives the results of some calculations of the reflection coefficient 
for electrons incident normally on a plane face of a metallic crystal. The physical 
situation is treated as being one-dimensional; and it is assumed that the potential 
energy of an electron is a sinusoidal function of distance inside the crystal, and 
obeys the classical image force law outside the crystal. The reflection coefficient 
is computed as a function of the energy of the incident electrons, over the range 
from to 20 electron volts, for a variety of values of the parameters which define 
the model of the crystal. 

1. Foreword 

THE work which is presented in this paper was undertaken as a result 
of conversations had with Dr. C. J. Davisson at various times during 
the years 1938 and 1939, when he was investigating the reflection of elec- 
trons impinging on the surface of a metallic crystal. The results for a simple 
special case of the general problem were published in 1939 1 . Thereafter the 
work on the general problem continued intermittently, and it was almost 
completed by the early part of 1942, when it was brought to a halt by the 
onset of wartime activities. Since then nothing has been done on the prob- 
lem, and the results already obtained have never been published in extenso. 
However, C. Herring and M. H. Nichols have included an illuminating dis- 
cussion of some of the more significant of the results in their recent mono- 
graph on thermionic emission 2 . 

Although the intervening years, by bringing new problems in physics to 
the fore, have caused this work to lose some of the interest which it possessed 
at the time it was being done, it still seems to be worth while to put the full 
results upon record. The present occasion, when his friends and former col- 
leagues are celebrating Dr. Davisson's seventieth birthday, is an especially 
appropriate one for this purpose. 

2. Formulation of the Problem 

We consider electrons moving with energy E and impinging on a plane 
face of a metallic crystal. (Fig. 1.) According to quantum mechanics there 
is certain probability R, generally neither nor 1, that an electron will be 
reflected by the crystal, and caused to move backwards toward the source; 
and there is the complementary probability 1 — R that the electron will 

1 Physical Review, v. 56, pp. 699-702. This paper will be referred to henceforth as 
[LAM, 1939]. 

1 Reviews of Modern Physics, v. 21, pp. 185-270 (1949). 
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penetrate the crystal, and flow away through the remainder of the circuit. 
We call R the reflection coefficient, and we can define it alternatively as the 
ratio of the intensity of the reflected electron beam to the intensity of the 
incident beam. 

We wish to calculate R as a function of E. In order to be able to do this 
effectively, it is necessary to idealize the actual physical situation quite 
drastically. (However, the idealization which we shall use preserves what 
seem to be the most important features of the physical situation.) On the 
other hand, once the idealization has been set up, the mathematical calcula- 
tions themselves will be carried through without approximations 3 . Hence, 




INCIDENT ELECTRONS 



REFLECTED ELECTRONS 



Fig. 1 — Reflection of an electron beam by a crystal. (Schematic representation). 



any discrepancies between the theoretical results and the results of experi- 
ment are to be attributed to the inadequacy of the model, and not to il- 
legitimate steps in the mathematical work. 

Our idealization of the physical situation can be described in the form of 
the three following assumptions: 

Assumption I. The problem may be treated as one concerning one-dimen- 
sional motion of electrons. Thus, we set up a rectangular coordinate system 
in space; and we assume that the crystal occupies the half-space X < 0, 
and that all of the point functions with which we are concerned depend 
solely upon the coordinate x. 

Assumption II. There exists a function V(x), such that an electron at the 
point x has potential energy V(x); and the behavior of an electron is gov- 
erned by the Schrodinger wave equation 

' Except, of course, simple arithmetical approximations, such as are involved in almost 
all calculations. 
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^ + k 2 [E - V(x)]f = 0. 



(1) 



(Here k 2 = Si^m/h 2 , where h is Planck's constant, and m is the mass of an 
electron.) This assumption deals in a summary way with various compli- 
cated processes involving electrons in crystals. Discussions of the validity 
of the assumption are to be found in various works on the electron theory 
of metals. 



S> a. 



v(x) 
o x 



Fig. 2 — Assumed potential energy as a function of the coordinate x. 



Assumption III. Specifically, the function V(x) is given by the formulae 

(2) 



V(x) = - V + Ki sin a(x - x ), 
x = e 2 /(4Fo), 



X < x 

X > x 



where e is the absolute value of the electronic charge, and V , V\ and a 
are suitable non-negative constants. (A graph of this function V(x) is 
shown in Fig. 2.) According to this assumption, an electron in the region 
x > .To is subjected to the classical image force. This is known to be in good 
agreement with the facts, at least if x is not too small.* Also, according to 
the assumption, the potential energy of an electron in the depths of the 
crystal is a periodic point function with a negative mean value. This part 
of the assumption is as correct as any assumption can be which attempts to 
account for the complicated actual processes in terms of a potential energy 
function. However, our particular choice of a periodic function is based 
largely upon mere considerations of mathematical convenience. Finally, we 

* See Herring and Nichols, footnote 2, p. 245 et seq. 
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observe that our V(x) is continuous, as physical considc...tions indicate 
that it should be. 

We can now state the mathematical problem before us in the following 
terms: 

V(x) being defined by (2), we are to obtain a solution \J/(x) of (1) satis- 
fying the following conditions: 

(a) In the region x > x the function i[f{x) exp (—2iriEt/h) represents an 
incident beam of electrons moving toward the left, and a reflected beam 
of electrons moving toward the right. 

(b) In the region x < Xq the average electron flow, if it is not zero, is direc- 
ted toward the left. 

(c) The function \j/(x) and its derivative \l/'(x) are everywhere continuous. 
Having obtained such a solution \/f(x) s we are then to compute the ratio of 
the intensity of the reflected electron beam to the intensity of the incident 
beam. In particular, we are to study the dependence of this ratio upon the 
quantities E, Vo, and Vy. 

The paper [LAM, 1939] already referred to dealt with the special case in 
which Vi = 0, i.e. the case in which V(x) is assumed to be constant in the 
region x < Xo. Consequently, we are now concerned chiefly with the cases 
in which Vi > 0. 

3. Generalities Concerning the Calculation of R 
In the region x > Xo the wave equation (1) takes the form 



g + *D+g,-a 



The general solution of this equation is of the form 

4,{x) = AH*) + BUx), 

where A and B are arbitrary constants, and ^i(.v) and fc(x) are two particu- 
lar solutions which we choose so that the functions ^i(.v) exp (—2iriEl/li) 
and ^(.v) exp (—2wiEt/h) represent beams of electrons, of unit intensity, 
moving to the left and right, respectively. 

In the region x < Xo the wave equation takes the form 

^4 + k 2 [E + Vo - Vi sin a(x - *b)]tf = 0. (3) 

ax 2 

We are concerned with a solution of this equation of the form 

4,{x) = cux), 

where C is a constant, and \l/3(x) is a particular solution such that the func- 
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tion ^ 3 (^) exp (—2iriEt/h) represents a state in which the average flow of 
electrons in the crystal either vanishes or is directed toward the left. 

The actual forms of the functions $i(x), \p2(x), faix) will be discussed 
presently. 

Now the continuity of the functions \J/(x) and yp'{x) gives us the system 
of equations 

>4^i(*o) + B$ 2 (x ) = Gp 3 (x ) 

Afo'ixo) + BtfM = W(*o), 

from which we can calculate the ratio B/A in terms of the ^i(xo), ^/(*o). 
Our required reflection coefficient R is | B/A | 2 , and so we obtain the for- 
mula 



R = 






(4) 



It was shown in [LAM, 1939] that the functions fa(x) and faix) are given 
by the formulae 

M*) = WU&, M*) = ^-a.i(-$), 

where 

£ = 2ikx&'\ X = -tlW/CSfi"*), 

and the symbols W\,\(£), PF_ x ,}(-£) denote the usual functions occurring 
in the theory of the confluent hypergeometric functions 4 . The earlier work 
gives us all the information concerning \f/i(x) and $i(%) that we shall require. 
Hence, in order to calculate R, we have, in effect, only to identify a suitable 
solution ypz(x) of equation (3), and then to calculate ^'(soV^fao). 

4. The Solution of Equation (3) 
In order to facilitate the use of known results, it is convenient to write 

«(* - *o) = 2z - ' *£ (E + V ) = el, i£ Vx - -26,. 

2 a* or 

Then equation (3) takes the form 

t| + (So + 25x cos 2*)* = 0. (3') 

This is one of the canonical forms of Mathieu's differential equation, for 

* E. T. Whittaker and G. N. Watson, "Modern Analysis" (Chapter XVI), Cambridge 
Univ. Press, 4th Ed., 1927. 
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which an extensive theory exists. We shall recall a few of the chief facts 
brought out in this theory. f 

Unless the constants O and 0i satisfy some one of certain special rela- 
tions, the general solution of equation (3') is of the form 

4, = K ie >"f{z) + K 2 e-»'f(-z), 

where /i is a constant determined by O and 0i, /(z) is a function which is 
periodic with the period t, and the K's are constants of integration. 

In certain ranges of values of the 0's, the constant y. is real, and in other 
ranges it is pure imaginary. When \i is real we can obviously take it to be 
positive; and then, in order that tyiix) may be bounded in the range x < x , 
we must choose \f/ 3 (x) to be the function e" z f(z). When \x is pure imaginary, 
we can take it to be i \ n | ; and then, in order that fa(x) shall represent a 
state in which the flow of electrons is to the left in the crystal, we must 
choose fa(x) to be the function e~" z f(—z). 

When n is pure imaginary we have a non- vanishing flow of electrons to the 
left in the crystal. Consequently, the intensity of the reflected beam must 
be less than the intensity of the incident beam. Hence, under this condition 
we must have R < 1. On the other hand, when ju is real there is no average 
electron flow in the crystal. Consequently, under this condition the inten- 
sities of the incident and reflected beams must be equal, so that R = 1. 
These considerations point to the importance of discussing, first of all, the 
conditions under which /z is real or pure imaginary. 

Figure 3 shows a well known diagram, modified slightly to suit our present 
purposes 8 . Here O and 0i are taken to be rectangular coordinates of a point 
in a plane, and the plane is divided into regions of two kinds (shaded and 
unshaded) by a system of curves. If the point (0o, 0i) is in the interior of 
one of the shaded regions, the above y. is real; if the point is in the interior 
of one of the unshaded regions, /z is pure imaginary. (If (0o, 0i) lies exactly 
on the boundary of one of the regions, we have a somewhat more compli- 
cated situation, which we do not need to consider here.) This diagram en- 
ables us easily to determine, for any fixed values of V and V\, the ranges of 
values of E in which we have R = 1 . We shall call these ranges of values 
of E the diffraction bands. 

Now our problem has been reduced to that of computing R for values of E 
which do not lie in diffraction bands. In treating this phase of the subject 
we shall follow the course of the actual calculations, without any examina- 
tion of ways in which the work might have been done more efficiently. 

t See, for instance, E. T. Whittaker and G. N. Watson, footnote 4, Chapter XIX. 
6 See, for instance, N. W. McLachlan, "Theory and Application of Mathieu Func- 
tions" (p. 40), Oxford University Press, 1947. 
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Of the many methods which have been devised for finding solutions of 
Mathieu's differential equation, the one which is conceptually simplest is 
that due to Bruns. This method can be described as follows: 




-2 -i 



0123456789 

Fig. 3 — Stability diagram for Mathieu's differential equation. 



10 11 12 13, 



Under the transformation 



\f/ = exp / <p(z)dz 



the Mathieu equation (3') goes over into the Riccati equation 

-f + <p 2 + 0o + 2ft cos 2z = 0. 
dz 



(5) 



We seek a solution of this equation in the form of a power series in the 
parameter 6\, say 

<p(z) = <p (z) 4- d!(pi(z) 4- 0i ^2(2) + ... , 

and we easily find that the functions <p (z), <Pi(z), ^2(2), ••• must satisfy the 
differential equations 



<po + <Po + 0o = 0, 

<pi + 2(poipi + 2 cos 2z = 0, 

<pz -\- 2<p<xp2 4" <pi = 0, 

<p 3 + 2<po<p3 + 2<W2 = 0, 

<Pi 4- 2<ptfpi + 2<pi<p3 + <P2 = 0, 



(6) 
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Obviously, in order that we shall get the solution which we require, <p{z) 
must reduce to — id Q when V\ = 0. Also p(z) must be periodic with the 
period 7r. These conditions, together with the equations (6), suffice to deter- 
mine <po(z), <pi(s), <p2(z), ■ • ■ successively and uniquely. We easily obtain the 
results: 

<po (s) = — id , 



2u -2 i 



(z) = - - r e - — 

** W " 2/ Ll — <9o 1 + 0oJ ' 



, v 1 r e iz 2 e~ iiz 



0o) 2 (2 - Go) (1 - 9t)9o (1 + 0„) 2 (2 + do) J 
M 1 r ^ , (4 - 3e )e 2i * 



+ 



* 4(2) = iMa~= 
+ 



16/ L(i - 0o) 3 (2 - 0„)(3 - e ) (i - 0„) 3 (i + e )e Q (2 - e ) 

(4 + 30„)<r 2 '' _ e - 6 " 

(i + e ) 3 (i - e )e (2 + e ) (i + e y(2 + O )(3 + a 
(n - 50 o y ,s 



j 



0o) 4 (2 - 0o) 2 (3 - 0o)(4 - O ) 
4(15 - 180o + 50oV'* 2(8 - 350 o 2 + 150 o 4 ) 



(1 - O ) 4 (1 + 0o)(2 - 0n) 2 (3 - o )0o (1 - O 2 ) 3 (4 - 2 )0, 



+ 



4(15 + 180o + 50b> 



2\ -4«2 



(1 + O ) 4 (1 - 9 )(2 + O ) 2 (3 + o )0o 

(11 + 50 o )g~ 8 " 

(1 4- O ) 4 (2 + 0c) 2 (3 + 0o)(4 + O ) 



] 



The functions <p b (z) and ^(z) have been computed also; but, because of 
their complexity, they will not be exhibited here. 

It is easily found that the expression ipai^/tiixo) appearing in equation 
(4) has the value <p(t/4)(ci/2) in terms of our present notation 6 . 

In principle we now have all the information we need to calculate the re- 
flection coefficient R- Furthermore, our experience showed that it is quite 
easy to compute R by this method, provided that the value of E does not 
lie too near an edge of a diffraction band. However, Brun's method proved 
to be unsuitable for calculating values of R for values of E in the neighbor- 
hood of a diffraction band edge, and we were forced to seek another method 
to obtain these values. After some tentative work with other methods, we 

6 VV'e must take account of the fact that the symbols i^(.v) and ip(z), which we are using 
for convenience, do not represent merely the same function with different arguments. 
In fact, after the change of the independent variable from x to z we have the following 
relation between the old and new xp's: [\p{x)} \ d = [>K*o + 2s/a — 7r/2a)] id = [<p{z)] aevl . 
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settled upon a method due to Whittaker, and this was used to complete the 
calculations. 

Whittaker's method is described briefly in Whittaker and Watson's 
"Modern Analysis," 4th Edition, p. 424. The method is developed more 
fully in papers by Ince 7 . We shall confine ourselves here to some summary 
indications of the nature of the method. 

The method leads to representations of the solutions of Mathieu's equa- 
tion by formulae which differ in structure depending upon the part of the 
(dl, 61) plane in which we are working. We shall give the formulae suitable 
for use in the neighborhood of the point 0o = 1, 0i = 0; the formulae for use 
in other parts of the plane are given by Ince. 

Given the values of 0o and 0i, we first determine a number a by means of 
the implicit equation 

2 6* 

e 2 = 1 4- 0! cos 2a 4- - 1 (-2 -h cos 4<r) - ± cos 2a 

Then we seek a solution of equation (3') in the form 
00 
xp = e" M 2 Kn+i cos[(2» + \)z - a] + b 2n+l sin [(2n + l)z - a]}, 

where n, the a's, and the b's are constants. We substitute the expression for 
\p into the differential equation, and determine values of the constants by 
imposing the condition that the resulting relation shall be an identity in z. 
After some rather intricate algebraic manipulations we finally arrive at the 
following results: 

0i • „ 30? . „ 30i . . , 

M = -sin2.--s 1 n2,-— sm4.+ ... 

ai = 0, 61 = 1 

a 3 - ^ sin 2a 4- M sin 4<x 4- | (- 2 -^ sin 2a 4- 9 sin 6cr) + • • • 

140? . „ . 440? . . , 
fl6= (8^ Sm2ff + (27K8^) Sm4a+ ''' 

350? . _ 

° 7= (iTJ8K^ Sin2 * + '" 

7 Monthly Notices, Royal Astronomical Society of London: v. 75, pp. 436-448; v. 76, 
pp. 431-442. 
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a. = 0(tf) 

fc-r+2«»>'+l(-7.t So " 4 ') 

+ §*( ~ 9 C ° S 2a + 7 C ° S 6a ) + " ' ' 

el ,40? _ . 6i /82 , 155\ . 

h = 77TS + 77^77^ COS 2(T + 7^7^ ( — COS 4(7 - — ) + 



192 ' (9)(8 3 ) (3)(8 4 ) \9 

(18)(8 3 ) T (12)(80 



*7 = , 4n * + TT^Ws COS 2ff + 



9 (180) (8 4 ) 



The calculated terms which are exhibited here enable us to calculate the 
solution \f/(z) to a certain accuracy, and this accuracy proved to be sufficient 
for our purposes. 

Although this method is very complicated analytically, it was found to 
be quite convenient for purposes of numerical calculation. 

5. The Reflection Coefficient for Large Values of E 

This work is concerned chiefly with the reflection coefficient for small 
values of E (actually up to 20 electron volts). However, it is interesting that 
we can obtain a simple approximate formula for R for indefinitely large 
values of E in the intervals between the diffraction bands. 

For this purpose we go back to Brun's method, and we write the dependent 
variable in equation (5) in the form ip = — ido + o. We find that the new 
dependent variable w satisfies the equation 

^ - 2*0 o w + w 2 + 20i cos 2z = 0, 
dz 

and we seek a solution of this equation in the form 

oji(z) w 2 (z) 



0) '= Wo 



(z) + -r^ + -5 s - + 



The functions oj„(z) are easily computed, and we finally arrive, in an entirely 
straight-forward way, at the result 

a a 

<*0r/4) = -fff, + -;+-!+•••• (7) 

0o fy 



- 
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In [LAM, 1939] we derived an approximate formula for R when E is 
large for the case in which V\ = 0. The work involved in that derivation, 
together with the equation (7), enables us to obtain the desired formula by 
a simple calculation. The result obtained is the following: 

^^L 1- ^ 1 .]- (8) 

The range of validity of the approximate formula (8) has not been de- 
termined. The nature of the derivation, and also the form of the result, leads 
us to suspect that the approximation is good only so long as the ratio 
V\/Vo does not exceed some bound depending upon the other quantities 
entering into the expression for R. The approximation certainly breaks 
down when Vi/V reaches the value 4Fo/(ae 2 ). However, this value is well 
above any of the values with which we deal with in this work. Consequently, 
we suspect that the formula can be used, to extrapolate our calculations of R 
to higher values of E, without serious danger of error in the cases which we 
consider here. 

6. The Calculated Results 

The reflection coefficient depends upon the independent variable E, and 
upon the three parameters V , V h and a. The effects upon R of taking vari- 
ous values of V and Vi seemed to be of greater interest than the effect of 
taking various values of a; and, consequently, we confined ourselves in the 
calculations to a single value of a, namely, a = ir X 10~ 8 cm -1 . This value 
of a makes the period of V(x) in the crystal equal to 2 X 10~ 8 cm. 

We took six values of V , proceeding in equal steps from 10 electron volts 
to 20 electron volts inclusive. These values adequately cover the range 
which is of interest in connection with actual metals. 

Including the calculations reported in [LAM, 1939], we have taken, for 
each of the values of V , five values of Vi, proceeding in equal steps from 
to 0.4 V Q . Although it is somewhat difficult to say just what value of Vi is 
most appropriate to the case of a specific actual metal, it appears that these 
values cover the range of values of interest adequately. 

The results of the calculations are shown in a self-explanatory form by the 
curves given in Figs. 4 to 9 inclusive. For the sake of unity, we have included 
the results which were previously published in [LAM, 1939]. 

7. Physical Discussion of the Results 

The results do not call for much discussion, especially in view of the dis- 
cussion which Herring and Nichols have given in the paper already referred 
to. However, there are a few observations which should be made. 
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ON REFLECTION OF ELECTRONS BY METALLIC CRYSTALS 905 

For fixed values of Vo and V u the reflection coefficient tends to decrease 
with increasing E over most of the range between any two successive 
diffraction bands. For fixed values of Vo and E, E being in a range between 
two diffraction bands, R tends to decrease with increasing V\. This is a re- 
sult which was not anticipated when the work was begun. As was expected, 
we find a tendency for R to increase with Vo when E and Vi/V are held 
fixed. 

The most interesting feature of the results is the behavior of R in the neigh- 
borhoods of the edges of the diffraction bands. Unfortunately, the range of 
values of E considered is not great enough to reveal this behavior very com- 
pletely. (The failure to consider a greater range of values of E was the result 
of our reluctance to embark again on the difficult numerical computations 
relating to the W±\,i(±£) functions. Since the necessary computations 
had been performed earlier for values of E up to 20 electron volts, we de- 
cided, unfortunately as it now appears, to confine ourselves to this range.) 

The behavior of the curves near the edges of the diffraction bands which 
occur in the neighborhood of E = 19 electron volts (when V is 18 or 20 
electron volts) does not require much discussion. The reader will observe a 
small dip in the curves for V x = V /10 just below these diffraction bands. 
The accuracy of the computations is believed to be high enough that we are 
justified in taking this dip to be entirely genuine. 

When Vo is 10, 12, or 14 electron volts the behavior of the curves for 
values of E in the neighborhood of zero is rather complicated. Herring and 
Nichols consider this behavior to be one of the most significant features of 
the results, and they have given a full discussion of it from the physical 
point of view.* In view of the availability of their discussion, we may con- 
fine ourselves here to a few brief remarks. 

In some of the cases which we are referring to now there are diffraction 
bands extending from E = to certain positive values of E. (These dif- 
fraction bands are shown in a self-explanatory way in the figures.) When 
such a diffraction band exists the complicated behavior of R for small val- 
ues of £ is a result of the existence of the diffraction band, and it is com- 
parable with the behavior of R in the upper part of the range of values of 
E in the case in which Vo = 20 electron volts. 

In the cases in which we do not have diffraction bands extending upward 
from E = we have to explain the complicated behavior of R in somewhat 
different terms. 

Assuming that we have such a case, let us momentarily ignore the fact 
that the physically significant values of E are non-negative, and consider 
E as an unrestricted real parameter. Under this convention concerning E, 

* Herring and Nichols, footnote 2, pp. 248-249. 
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we find that there is a diffraction band lying in the range of negative values 
of E, and extending more or less near to the point E = 8 . We shall call this 
a fictitious diffraction band. Now it is immediately clear that the complicated 
behavior of our curves arises from the fact that the small positive (and hence 
physically significant) values of E concerned are near the upper edge of the 
fictitious diffraction band. 

8 Specifically, what is meant is that there is such a range of values of E in which the 
exponent n is real. 



